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THE MOBPHOLOGY OF ENERGETIC o"^ IONS 
DURING TWO MAGNETIC STORMS: TEMPORAL VARIATIONS 


R. D. Sharp, R. G. Johnson, and E. G. Shelley 
Lockheed Palo Alto Research Laboratory 
Palo Alto, California 9^304 


ABSTRACT 

• 4 " 

A study has been conducted of the morphology of precipitating 0 and 
H ions in the energy range 0.7 ^ E 12 keV during the storm- time period 
from December l6-l8, 1971, which encompassed two principal magnetic storms. 
This paper describes some of the results of this study with emphasis on 
the temporal variations of parameters characterizing the intensity, average 
energy, and spatial location of the zones of precipitation of the two ionic 
species. One of the principal results was the finding that the intensity of 
the precipitating 0^ ions was well correlated with the geomagnetic indices 
which measure the strength of magneto spheric substorm activity and the 

■ 4 * 

strength of the storm-time ring current. Since the 0 ions are almost 
certainly of ionospheric origin the correlations indicate that a previously 
unknown strong coupling mechanism existed between the magnetosphere and the 
ionosphere during the storm period. Some other morphological features appar- 
ent in the data cure: l) the storm-associated initial increase of the 0^ 

ions on the nightside (0300 LT) was found to lead that on the dayside (1500 
LT) and lag the initial nightside increase by more than une hour in both 
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storms; 2) a strong correlation was observed between the venriations in loca- 
tions of the 0 and H precipitation regions on both the day and nightside 
crossings; and 3) the average energies of the 0^ and precipitations were 
significantly correlated on the dayside crossings. The implications of 
these results with respect to some phenomenological models of the 0^ morph- 
ology are discussed. 

The total worldwide precipitated ion energy flux has been estimated 
during the period of the study and compared to the ring current energy 
content as measured by Dst. The comparison indicates that ion precipita- 
tion was an important loss mechanism for ring current energy during the 
December 17-18, 1971 magnetic storm. 


INTRODUCTION 


The discovery of 0^ ions with energies of up to 12 keV in the magneto- 
sphere was reported by Shelley et al . [1972]. Synoptic studies of several 
aspects of the morphology of these ions were described by Sharp et al . [197^], 
Shelley et al . [197^] » and Johnson et al . [1975]- The energetic 0^ ions were 
inferred to be of ionospheric origin. They were observed in every storm 
studied over a one-year period and at all local times. They were also ob- 
served at reduced intensities during non- storm times. This is the first of 
two papers which will present more detailed results from a statistical study 
of the temporal and latitudinal variations of the ions emd the accompany- 
ing protons in the same energy range (0.7 ^ E < 12 keV) during two magnetic 
storms. This paper will present integral parameters, each characterizing 


- 4 - 


a single traversal of the precipitation zone. Paper II will describe the 
latitudinal variations averaged over each of the two storms. A simileur 
statistical analysis of a comparable body of data on auroral electrons 
€Uid the more energetic component of the auroral protons has been pub- 
lished [ sharp and Johnson , 19^8; Sharp et a' . , 1969]. 

The two storms studied here occurred in the period December l6-l8, 
1971. Both were initiated with sudden commencements which occurred at 
1904 UT on December 16 and l4l8 UT on December 17. The first storm had 
only a small main phase with a peak Dst of 54v. The second storm had a 
classic main phase with a peak Dst of 17IY and a recovery phase lasting 
until about 2300 UT on the l8th. The relevant geomagnetic indices charac- 
terizing this period ?re shown in Figure 1. Extensive data on the equa- 
torial particle distri jutions during this same period were acquired by 
the Explorer 45 satellite [ Smith and Hoffman , 1973; Hoffman , 1973; Williait)R 
and Lyons , 1974a, b; Lyons and Williams , 1975; and Wi lliams , 1975]. 

The data presented here were obtaintd from an energetic ion mass 
spectrometer experiment on the low-altitude polar satellite 1971-089A. 

The experiment contained three spectrometers. Each consisted of an elec- 
trostatic analyzer in series with a crossed field velocity filter which 
provided both mass- per -unit-charge and energy-per-unit-charge information 
on the measured particles. The satellite was in an approximately circular 
orbit of 800 km altitude and 93° inclination. It was in the local time 
plane 03OO-I5OO during the period described here. The satellite was sta- 
bilized about three axes and the spectrometers were oriented at 55° to 
the local zenith so that they generally sampled precipitating particles 
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during their traversals of the auroral and sub-auroral regions. The nine 
energies per unit charge sampled by the experiment were 0.74, 1.01, l,4l, 
2.14, 2,92, 4.07, 6.3, 8.6, and 12.1 keV. Mass-per-unlt-charge scans at 
each of these values were obtained In 6.l4 seconds. Representative mass- 
per-unlt charge distributions are given In Sh elley et aL [1972] which also 
contains a more detailed description of the experiment and preliminary 
results from this same storm-time period. 

ANALYSIS AND RESULTS 


In order to investigate the variations in the properties of the ion 
fluxes on a time scale appropriate to the period of the storm, we have 
formulated a group of parameters based on the average or integral fluxes 
during a single complete traversal of the northern precipitation zone. 

This gives us a measurement in each of the northern precipitation zonv.6 
(day and night) about every 100 minutes which is comparable to the time 
resolution of the various geomagnetic indices characterizing the disturbed 
period. To formulate these parameters for both the 0^ auid ions we have 
summed the counts in the region of peak response in the mass-per-unlt-charge 
sweep for each ion species and subtracted the background from an equivalent 
number of channels. Then for each six-second measurement cycle, we have 
computed the integral energy flux and average energy for each species in 
the energy range 0.7 ^ E s 12 keV. For the intensity parameter, we have 
integrated the flux over the range in invariant latitude 44° ^ s 8o°. 

For the hardness parameter we have computed the average energy of each ion 
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specles over the same interval. To characterize the location of the pre- 
cipitation zone we have computed the invariant latitudes at which 10^, 50^ > 
and 90^ of the energy flux latitudinal Integrals were reached. 

The results for the northern hemisphere observations during the period 
1200 UT on De'iember l6 to 2400 UT on December 18, 1971 are presented in 
Figures 2 through 5* Figure 2 shows the integral energy flux values. They 
he e been multiplied by the satellite velocity and corrected for the (gener- 
ally small) deviations of the trajectory from meridians of magnetic longi- 
tude so uhat they are an approximation to the Insteuitaneous value of the 
energy flux precipitating through a one-cm-wlde slit aligned perpendicular 
to the precipitation zone. The error bars represent counting statistics 
only. There are additional experimental uncertainties, discussed by Shelley 
et al . [ 1972 , 1975 ], which we estimate as approximately + 

One sees in Figure 2 that there is a general correspondence between 
the 0^ zone integrals and the magnetic activity indices shown in Figure 1 
with large flux increases beint observed shortly after the sudden coaaaence- 
ments. The nigbtside (0300 LT) 0^ increase leads the dayside ( 15 OO LT) in- 
crease and lags the increase by over an hour in each storm. On the night- 
side, the 0^ and intensity levels are generally comparable while on the 
dayside the 0^ precipitation is considerably weaker than the precipitation 
in this energy range. 

Figures 3 snd 4 show the zone average energy values. For Figure 3, we 
defined the average ion energy as the ratio of the zone integral energy flux 
divided by the zone integral number flux. For those cases where the latter 
was less than two standard deviations above zero the average energy values 
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were not considered significant and were deleted. The error bars represent 
counting statistics only. This definition for the avere^e ion energy tends 
to give the tnost weight to the intense precipitation events in each zone 
crossing. In order to evaluate the significance of the observed trends we 
defined the alternative hardness parameter shown in Figure h. fbr this 
parameter we computed the average energy for each Individual (six-second) 
spectral measurement that was considered statistically significant and then 
computed the grand average value of these quantities over the precipitation 
zone, weighing each measurement equally. The significance criterion chosen 

was that the integral number flux in the individual measurement had to be 

% 

greater than twice its standard deviation. The error bars in this figure 
represent the standard deviations of the means of the.'«e individual average 
energies and not counting statistics. For those cases where less than ten 
individual measurements were included, the number of measurements has been 
indicated on the figure. In examining Figures 3 and U we see that the aver- 
age energies do not vary substantial ly between the two ion species ov be- 
tween the day and nightside measurements. There is a correspondence in the 
temporal variation!.' of che 0^ and values on the dayside that is apparent 
in both Figures 3 and U so that we feel it is not the result of biasing due 
to a specific weighing technique. We will discuss the implications of this 
observation in the next section. 

Figure 5 chows the locations of the precipitation zones. The center of 
the region, measured by the Invariant latitude of the 50^ point of the zone 
integral energy flux, is indicated with a circle for 0^ ions and a square 
for ions. The horizontal bars represent the extent of the region as 




FIGURE 5 
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measured by the 10^ and 90^6 points in the same parameter. In several in- 
stances during this period the proton precipitation extended above 80° 
invariant latitude. The satellite orbit on some passes did not extend 
higher than 80° so in order to form a uniform data base we truncated the 
zone integrals at that latitude. Thus for those points in Figure 5 in 
which the upper bar is in the vicinity of 80°, that point does not neces- 
sarily mark the upper latitude of the precipitation region. Despite this 
distortion we see a detailed tracking of the precipitation zones of the two 
types of ions on both the day and nightsides. There is a latitudinal dis- 
placement between the two species with the 0^ zone generally located equat- 
orward of the zone by several degrees. There is a correlation between 
the results at the two local times as well as between the two ion species. 

DISCUSSION 


As noted above in connection with Figure 2 the observed nightside 0^ 
fluxes are more intense than those on the dayside and there is a time delay 
in both storms between the initial 0 increases on the night and daysides. 
This is su gg estive of a longitudinal drift model with the principal injec- 
tion occurring on the nightside. In the energy range of these measurements 
longitudinal drift is strongly a.fected by the corotation and convection 
electric fields. Shelley et al . [1974] showed from a study of the local 
time dependence of the 0^ ions that the median peak precipitated flux is about 
an order of magnitude higher in the LT interval 21-09 hours than it is in 
the interval 09-21 hours. The asymmetric location of the region of high 
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median intensities vlth respect to the midnight meridian is perhaps related 
to longitudinal drift since, according to most magnetos pheric electric field 
models, there is a net eastward drift expected for trapped ions with energies 
of a few keV or less in the dawn section in this L range [ Kavanagh et al . , 
1968; Mcllwain, 1974]. If one looks at this hyEK,;.hesis more queuititatively, 
however, one finds that the observed delay times are probably too short to 
be explained by a longitudinal drift of the ions from the local time 
sector of high-median peak intensities (21-09 hours). Consider first the 
D-^^cember 17 storm for which the more definitive early-time dayside data 
were obtained. The time delay between the initial increases observed on 
the night and daysides is 1.4 hours. Depending upon the details of the 
magnetospheric electric field model and the specific particle energy, ions 
in the few-keV range injected in the midnight sector can longitudinally 
drift either eastward or westward. The shortest delay time expected for 
eastward drifting ions would be for those of 0 enargy and this would be 
about five to ten hours for drift from 0900 to the observation point at 
1500 LT [ Kavanagh et al ., I968; Mcllwain , 1974]. The shortest delay ex- 
pected for westward drifting ions would be for the most energetic particles. 
Assuming E = 7.8 keV (the average energy of the Initially observed dayside 
group) and L = 3.8 (the position of the median integral energy flux for this 
group) and neglecting convection for a lower limit we find that about nine 
hours would be required to drift from 2100 to 1500 LT for ions mirroring 
at low altitudes. For equatorially trapped ions, the corresponding time 
is about 6.5 hours. For the December I6 storm, there are possible ambi- 
guities because of the low flux intensities and data gaps, but the time 


- 15 - 


delay between the Initial increase on the night and dayaides is most prob- 
ably « 3.2 hours. A similar d^lft calculation to that described above shows 
that a delay time of greater than 6,3 hours would be required for ions mirror- 
ing at low altitudes and 4.U hours for equatorially mirroring ions. We con- 
clude that the initially observed 0^ ions in both storms are not likely to 
be convected around from the nightside but more probably are accelerated 
locally. The observed delay times can perhaps be understood in terms of the 
time required for the energetic ring-current protons (a proposed energy 
source for the 0^ acceleration [ Cladis , 1973a, b]) to drift around westward 
from the nightside sector. For example, equatorially- trapped 50-keV protons 
would drift from midnight to 1500 LT at L = 6 in about one hour. 

If the initially observed 0^ ions on the dayside are indeed accelerated 
locally, this provides a significant constraint on models of the acceleration 
mechanism. It must be one that operates effectively under widely different 
ionospheric conditions to accelerate particles from both the day and night- 
side ionospheres. 

Although the initially observed 0^ ions have apparently not convected 
from the nightside, those observed later in the storm might have done so. 

The observed softening in the average energy of the dayside 0^ ions about 
five hours after the initially observed nightside increase in the December 
17 storm is in about the right period for the expected time of arrival of 
eastward convecting soft ions. 

We have remarked on the general correspondence between the 0^ intensity 
parameters and the geomagnetic disturbance indices. To exeunine this rela- 
tionship in more detail we have superimposed plots of the Dst and AE indices 
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and the nightside 0^ intensity parameter on the same logarithmic scale in 
Figures 6 and 7. For the Dst plot (Figure 6) we have inverted the ordinate 
so that increasing ring current intensity corresponds to increasing 0^ intens- 
ity. One sees that there is indeed a significant correlation between the 
various quantities throughout the period of the study indicating that the 
energized 0^ results from a strong coupling between the ionospheric source 
and the magnetosphere. The correlation with Dst is suggestive of a ring- 
current energy source as has been suggested by Cladis [1973a, b] and Brice 
and Lucas [1975], or alternatively that the 0^ ions contribute substantially 
to the ring current. On the other hand the correlation with AE suggests the 
possibility of a source associated with magneto spheric substorms, for example 
through the acceleration of the ionospheric ions by the field aligned elec- 
tric fields expected to result from the anomalous resistivity caused by the 
enhanced Birkeland currents associated with substorms [ Thorne , 1975] • ions 
have been reported during substorms by Johnson eo al . [1975] • If the 0^ sourct 
mechanism were indeed proportional to the ring-current intensity (Dst) or the 
strength of electrojet activity (AE) and the 0 ions were directly precipi- 
tated with negligible trapping lifetime, we would expect a detailed corres- 
pondence between the precipitated 0^ intensity and the appropriate geomagnetic 
index. We see in Figure 6 that the ring-current hypothesis fits reasonably 
well to the early-time data in both storms. The apparently shorter decay time 
for the 0^ intensity than for Dst and the large increase in 0^ intensity near 
1500 UT on December l8 which is not reflected in Dst argues against this model 
for the later phases of the storm. 
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The substorm hypothesis relevant to Figure 7 also gives a rough quali- 
tative fit to the observations. In particular, we note the corresponding 
peaks in the two paraaeters around midday on the l8th. We would, however, 
hesitate to conclude that the correlation with AE shown in Figure 7 is defin- 
itely superior to that with Dst in Figure 6 on the basis of this limited 
amount of data. 

A model utilizing a ring-current energy source that would allow for 
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different decay times for the observed 0 intensity and the Dst index is 
one in which the observed 0^ is precipitating from a trapped parent popula- 
tion which is controlled by decay processes different from those for the 
energetic protons which presumably are primarily responsible for the ring 
current, at least at later times in the storm [ Berko et al . , 1975]. 

The best estimate of the 0^ decay time to use for these considerations 
is probably obtained from the data in the OOOO-O60O UT time period on December 
18 when both AE and Dst are declining nearly monotonically and injection pro- 


u eases are presumably weak. 
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precipitation is falling off with about a 2-1/2 hour time constant. This is 
substantially shorter than the estimated decay time resulting from strong 
pitch-angle diffusion or charge exchange for a longitudinally uniform trapped 
population of 0^ ions. For the nightside 0^ data, the invariant latitude of 
the median integral energy flux varies from 64° to 70° during this period 
and the average 0^ energy is about 4.5 keV. The minimum lifetime for strong 
pitch-angle diffusion corresponding to these values for 0^ ions varies over 
the range from approximately 17 to 78 hours [ Kennel , I969] . The appropriate 
charge exchange lifetime for equator ially mirroring 0^ ions is greater than 
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four days baaed on the lifetimes for proton charge exchange given by Swisher 
and Frank [1968] modified by the ratio of the charge exchange cross sections 
for H on H and 0 on H summarized in Fite et al . \\So2'\. 

Let us consider two simple trapping models that attempt to explain 
both the observed 0^ decay time and the asymmetry in the local time distri- 
bution of peak 0^ intensities reported by Shelley et al . [1974], assuming 
a symmetric injection region around midnight. For the first model we assume 
that the 0 ions are trapped but precipitating under strong pitch-angle diffu- 
sion so that the time dc?pendence of the 0^ intensity parameter also repre- 
sents the time dependence of the trapped flux intensity at 0300 LT. Consider 
the decay of a longitudinally limited cloud of ions injected into the 2100- 
0300 LT sector. The lifetime of the population at 0300 LT is determined by 
both precipitation and the depletion of the trapped population by longitudinal 
drift and convection. A detailed convection calculation is beyond the scope 
of this work but we can get some idea of the time scales involved by consider- 
ing the motion of zero energy ions using current magnetospheric electric field 
models [ Mcllwain , 1972, 1974], A cloud of such ions with the local time extent 
assumed above would take of the order of an hour to convect past the observa- 
tion point at 0300 LT. More energetic ions would take longer because their 
gradient and curvature drifts oppose the convection motion. Depending on 
the details of the assumed electric fields, ion energies, and pitch-angle 
distributions and including the effect of precipitation, the model could prob- 
ably be fit to the observed 0^ decay time. 

The decrease in the local time distribution of peak flux intensities at 


sa 0900 LT reported by Shelley et al . [1974] might be interpreted in this 
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model aa resulting from the decay of the trapped population through precip- 
itation in the time it takes to convect from the injection region to O 9 OO LT. 
For this to occur the convection time would have to be comparable to the 
precipitation lifetimes and this seems somewhat long, but not clearly ex- 
cluded. Some local acceleration on the dayside would still be required in 
this model to explain, for example, the small time delay observed between 
the rise of the dayside and night side 0^ fluxes shown in Figure 2. 

For our second model, let us assume again that the observed 0^ ions 
are precipitating from a parent trapped population but that the precipita- 
tion is not necessarily proportional to the trapped flux intensity. We 
ceui assume, for example, that during the recovery phase the strength of 
pitch-angle diffusion declines and the observed 2-1/2 hour time constant is 
not representative of the decay of the trapped flux intensity but of the 
precipitation mechanism. The extent toward morning of the local time asym- 
metry could then be determined by how far the 0^ ions convect in the time 
period when the precipitation mechanism is most effective. In the storms 
described here, the time from the initial 0^ increase until the flux begins 
to fall off more rapidly than Dst is about 10 to 12 hours. Again, this 
is getting into the range where a fit to a convection calculation appears 
possible depending on the details of the model. 

The correlations between the protons and the 0^ ions apparent in 
Figures 3 to 5 apply additional constraints to these phenomenological 
models. The tracking of the location of the two precipitation zones could 
arise from common motions of trapped populations of the two species; from 
a process in which the ions were energized by a mechanism involving the 
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trapped protons as an energy source; or from a process in which the protons 
and 0^ ions were accelerated together from the ionosphere by the same mech- 
anism. In the latter two cases a substantial trapped 0^ population would not 
necessarily be involved. The correlations in the variations in average 
energies of the two species on the dayside seem most directly interpret- 
able in terms of storm-associated acceleration processes acting upon the 
parent trapped populations, or with a common ionospheric source and a 
common acceleration mechanism varying in efficiency during the storra-time 
period. This latter interpretation would have far-reaching effects on the 
present state of thinking about the origin of the energetic magneto spheric 
proton population [ Hill , 1974]. 

We have briefly explored a few of the constraints which these morpholog- 
ical results can exert on possible source and energization mechanisms for the 
0^ ions. It is hoped that more quantitative models will be investigated by 
others utilizing the experimental results presented here and in Paper II. 

As discussed in the "Analysis and Results" section, the zone integral 
intensity values can be considered as an approximation to the energy flux 
precipitating into a one-cm-wide strip transverse to the precipitation zone. 

If we knew the local time distribution of the fluxes, we could estimate 
the instantaneous worldwide precipitated ion flux in this energy range dur- 
ing the course of the storm. The local time distribution of precipitating 

protons can be estimated from the observations of Bather and Mende [l971] 

P 

and Mende (private communication). A rough estimate for the local time dis- 
tribution of the 0^ ions can be obtained from the work of Shelley et al. 


[ 1974 ]. From these results we find that a reasonable approximation to the 
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local time average of the flux intensity can be formulated from the average 
of our measured values at 1500 and 03OO hours local time. By assuming con- 
Jugacy and Isotropic pitch-angle distributions over the upper hemisphere, 
we can approximar.e P, the total worldwide precipitated energy flux in the 
range (0.7 ^ E ^ 12 keV) as 

2 

P « 2tt Rg (Fjj cos + Fjj cos ^^) ergs/ sec 

where F = the zone integral energy flux (Figure 2), A = the inveuriant lati- 
tude of the 50j& point of the zone integral (Figure 5), and the subscripts 
D and N refer to day (1500 LT) and night ^300 LT) , respectively. Plots of 
P for both the and 0^ ions are given in Figure 8. 

If we assume that the precipitating ion fluxes originate from (or derive 
their energj' from) the trapped population responsible for the stormtime ring 
current we can ask whP^. fraction of the ring-current energy is lost through 
ion precipitation. ' j instantaneous value of Dst can be related to the 

total energy content ox' the trapped ring-current particles (E) by the expres- 
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Sion Dst(v) = 2.6 X 10 E, where E is given in ergs [ Sckopke , I966; Frank , 

19671. 

In order co have directly comparable quantities we would like to make 
the comparison between P and changes in Dst in a time period when ring-current 
Injection processes are negligible. Da-/ie and Parthasarathy [1967] have shown 
that during such periods AE is relatively low and Dst decays exponentially 
with a time constant of about ten hours. The time period from 0000 to 080O U1 
on December I8 fulfills these conditions. During this period 2.8 x 10 ergs 
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«re lost from the ring current according to the above relationship. By 
integrating the curves in Figure 8 over this interval, we find that the 
worldwide energy loss from precipitating ions with energy between 0.7 and 
12 keV is 6 X 10 ergs of 0 ions and 6 x 10 ergs of protons, thus 
accounting for about % of the lost ring-current energy. 

The major portion of the ring current energy is carried by protons 
with energies above the range oi our spectrometer. For a first-order 
estimate of the energy loss in this portion of the spectrum, we can 
utilize the Explorer 45 measurements of trapped protons in the energy 
range from 1-500 keV taken during this same time period at the latitude 
of the maximum of the equatorial energy density of the peurticles (L = 1 . 6 ). 
Integrating the orbit 102 spectrum in Figure 3 of Smith and Hoffman [1973] 
we find that only of the energy density is carried by protons with ener- 
gies less th ^n 12 keV. If we assume that this spectrum is applicable to 
both the precipitating and 0^ ions we conclude that about 6o^ of the 
ring-current energy loss can be accounted for by precipitating ions. If 
we assume that the published spectrum is applicable only to the precipitating 
tons and that the spectr ’m of the precipitating 0^ ions cuts off sharply 
at 12 keV, we still can account for about 403^ of the ring-current energy 
loss in this way. It should be noted that the Explorer 45 spectrum aver- 
aged over the energy range 1 i E ^ 12 keV gives an average energy which is 
consistent with the average energy of the ions measured during this period 
with the 1971-O89A experiment. 

The Explorer 45 energy spectrum which we have utilized was constructed 
with data from two types of instruments. The lower energy piotion (E < 30 
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kev) was measured with an electrostatic analyzer which would have counted 
any 0 ions in the equatorial plane. The upper energy portion was measured 
with solid-state detectors which are relatively insensitive to 0^ ions be- 
cause of the increased lossen of the heavier ions in the dead layer of the 
detector. By utilizing the spectrum as published we have essentially 
assi^ed that the equatorially trapped 0^ fluxes are negligible during this 
period. An alternative assumption is that the fluxes have a similar 
equatorial pitch-angle distribution to that of the protons. In this case 
we would modify the published spectrum at the low energies by deleting the 
0 contribution according to the H /O ratio measured at low altitudes 
during this period. 

Under this assumption only 6^ of the trapped proton energy density 
Is carried by protons with energies below 12 keV and either about 80^ or 
50 ^ of the ring current energy loss is accounted for by precipitating ions 
during this period according to the two assumptions on the 0^ spectrum dis- 
cussed above. 

As indicated, these are only meant to oe rough estinutes. It appears, 
however, that under the stated assumptions we can conclude that ion precip- 
itation was an important loss mechanism for the ring-current energy during 
this storm. If a large fraction of the precipitating ions derive their 
energy from sources other than the ring current, this ccnclusion is of 
course not valid. One obvious such source is the magnetosheath which is 
known to contribute to the dayside proton preo-i pitation at high latitudes 
[ Heikkila and Wlnningham , 1971]. We do not feel that the polar cleft 
particles contribute substantially to the measured integrals how-v-'^r; rinco. 
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as will be seen in Paper II, the average proton energies even at high lati- 
tudes on the dayside remain in the few-keV range, v?aa.l above the average 
proton energy in the cusp. In any case, even deleting the dayside proton 
data completely would only reduce the above estimates by about 25?^. 

Our conclusions about the importance of ion precipitation as a loss 
mechanism in this storm are in disagreement with the conclusions of Mizera 
[197U3 for the March 19^20, 1969 magnetic storm. The details of his cal- 
culation were not presented, but on the basis uf electrostatic analyzer 
data he concluded that less than 1^ of the ring-current energy lost during 
that storm could be attributed to precipitating protons with E > 12 keV. 

SUMMARY AMD CONCLUSIONS 

Some of the principal results of thlt. study of the zone integral 
parameters of the precipitating 0^ and H*' ions during the December 16-18, 
1971 storm period are: 

1) The intensity of the precipitating 0 ions was found to be well 
correlated with the geomagnetic indices which measure the strength of 
magneto spheric substorm activity and the strength of the storm-time ring 
current. Since the 0^ ions are almost certainly of ionospheric origin, 
these correlations indicate that a previously unknown strong coupling 
mechanism existed between the magnetosphere and the ionosphere during the 
period of the study. 

4. 

2) The storm-associated initial increase of the 0 ions on the 
nightside (0300 LT) was found to lead that on the dayside (1500 LT) and 
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lag the initial nightside increase by more than one hour in both storfus. 

A consideration of ionic transport processes leads to the conclusion that 
the taechanism accelerating the 0^ ions is probably operative in the dayside 
as well as the nightside ionosphere. This is significant in that it implies 
that this unknown mechanism can be operative over a wide range of ionospheric 
and magneto spheric conditions. 

3) Correlations have been found between the locations of the 0^ and 
precipitation zones and between the average energies of the two ionic 
species. From these and other morphological features, it is inferred that 
the 0^ ions are probably either accelerated together with the protons or 
coexist with them as a trapped population for an extended period. Either 
of these hypotheses imply that the ionospheric contribution to the storm- 
time ring current may be more significant than has previously been consid- 
ered. 

U) The total worldwide precipitated ion energy flux has been estimated 
during the period of the study and compared to the ring-current energy con- 
tent as measured by Dst. The comparison indicates that ion precipitation 
was an important loss mechanism for ring-current energy during the 17-18 
December 1971 magnetic storm. 
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